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Abstract. Using a recent path integral representation for the T-matrix in nonrelativistic potential scat- 
tering we investigate new variational approximations in this framework. By means of the Feynman- Jensen 
variational principle and the most general ansatz quadratic in the velocity variables - over which one has 
to integrate functionally - we obtain variational equations which contain classical elements (trajectories) 
as well as quantum- mechanical ones (wave spreading). We analyse these equations and solve them numer- 
ically by iteration, a procedure best suited at high energy. The first correction to the variational result 
arising from a cumulant expansion is also evaluated. Comparison is made with exact partial-wave results 
for scattering from a Gaussian potential and better agreement is found at large scattering angles where 
the standard eikonal-type approximations fail. 

PACS. XX.XX.XX No PACS code given 



1 Introduction 

Variational approaches to quantum mechanical scattering 
have a long history and are well covered in standard text 
books (e.g. [T1[5J|3J|3]). Actually, as formulated in ref. '5 , 
it is "... possible to construct systematically a variational 
principle for just about any given quantity of interest, pro- 
vided that the entities which enter into the definition of 
(that quantity) are uniquely defined by a set of equations 
For quantum physics the observables obtained from a 
solution of Schrodinger's equation are of particular inter- 
est and therefore variational principles have been available 
since the beginning of quantum mechanics. Best known 
is the Rayleigh-Ritz variational principle for the ground- 
state energy of bound systems but the continuous spec- 
trum is also accessible to a variational treatment. Most 
prominent among the variational principles for scattering 
are Schwinger's and Kohn's IJ^ expressions. In particu- 
lar, Kohn's variational principle is widely used in nuclear 
[51IMTU] and atomic jllyl2j physics for an approximate 
description of few-body scattering near thresholds. 

While these approaches benefit from the flexibility which 
ingenious trial wave functions offer it is well known that 
in many-body systems or in field theory the use of wave 
functions (or functionals) ceases to be useful. The path- 
integral method where one integrates functionally over the 
degrees of freedom weighted by the exponential of the clas- 
sical action is much more general although the cases where 
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one can actually perform the path integral are rare. There- 
fore, in general, one has to resort to approximations, such 
as perturbation theory or brute-force numerical evalua- 
tion of the functional integral on a (space-time-) lattice. 
If this is not appropriate or feasible, one may use a vari- 
ational principle extended to actions. A prime example is 
the Feynman- Jensen variational principle which has been 
used to obtain the best semi-analytic ground-state energy 
of an electron in an ionic crystal (the polaron problem 

Oddly enough, scattering has mostly remained outside 
the path-integral approach and it is only at zero energy 
that bounds for the scattering length have been obtained 
from the path integral in the imaginary-time formula- 
tion [14'T5 . Expanding on previous attempts [16' a real- 
time path-integral representation for the nonrelativistic 
T-matrix in potential scattering has recently been derived 
[17] . In this formulation the particle travels mainly along a 
simple reference path while quantum fluctuations around 
this path are taking into account by functional integration 
over velocities. It has been shown that this description 
gives the exact Born series to all orders if an expansion in 
powers of the potential is done and reduces to the eikonal 
approximation (valid at high energies and small scatter- 
ing angles) when the quantum fluctuations are neglected 
altogether. Taylor-expanding the action around the refer- 
ence path and performing the Gaussian functional inte- 
grations term by term, a variant of the systematic eikonal 
expansion of the scattering amplitude derived by Wallace 
[IH] is obtained (higher orders have been calculated in ref. 
[19)1. This is very promising for applications to many- 
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body scattering as the eikonal approximation is the basis 
of Glauber's very successful theory of high-energy scatter- 
ing from composite targets p0,21 . 

Given this affinity to a geometric description of high- 
energy scattering and the success of Feynman's treatment 
of the polaron problem it seems interesting to study how 
a variational approach to potential scattering performs in 
this framework. With a simple (linear) ansatz for the trial 
action this has been investigated in ref. |22) where it was 
found that the classical trajectory - and not a straight- 
line path as in the eikonal approximation - determines the 
scattering dynamics^. Numerically, very promising results 
in potential scattering have been obtained in cases where 
the eikonal expansion fails. 

It is the purpose of the present work to generalize this 
work by allowing for the most general quadratic + linear 
trial action. One may expect that the additional quadratic 
term describes the wave-spreading characteristic for the 
exact quantum theory thus leading to a much better de- 
scription of the scattering process. 

To be in agreement with the high-energy eikonal ex- 
pansion such an ansatz must allow for anisotropic terms 
which have already been shown to improve a variational 
calculation in a scalar field theory [21] . Since the Feynman- 
Jensen variational principle is the first term of a cumulant 
expansion it is also possible to calculate systematic cor- 
rections. We do it here by evaluating the second cumulant 
which is similar as in the polaron problem P5ll26j but also 



where 



more challenging as we have to deal with the complex 
scattering amplitude. 

The paper is organized as follows: In sect. 2 we present 
the essentials of the path-integral representations of the 
T-matrix in potential scattering so that we can apply the 
Feynman-Jansen variational principle in this setting. Sect. 
3 contains the variational ansatz and derives the ensuing 
variational equations. Some properties of the solutions and 
special cases are then discussed and the correction by the 
second cumulant is given. Sect. 4 then presents our nu- 
merical results for high-energy scattering from a Gaussian 
potential and comparison with exact partial- wave calcula- 
tion of the scattering amplitude as well with other approx- 
imations discussed in the literature. The work concludes 
with a summary and outlook for further work and appli- 
cation. Most of the technicalities as calculation of various 
path-integral averages and numerical details are collected 
in four appendices. 



2 Path-integral representations of the 
scattering amplitude 

Recently two variants of a path-integral representation for 
the T-matrix in potential scattering have been given |17) 
in the form 



= i 



K 



5(b) - 1 



(2.1) 



^ Corrections to the straight-line trajectory also have turned 
out to be important in heavy- ion collisions [23) . 



K = -(k,-Kk/) , |K| =A:cos- (2.2) 



q = k/ - ki , g EE |q| = 2fc sin - 



(2.3) 



are the mean momentum and momentum transfer, respec- 
tively. fc^/(2m) is the scattering energy (we set h — 1) 
and 9 the scattering angle. For lack of a better nomencla- 
ture we will call S{h) the "impact-parameter 5-matrix" 
although eq. (j2.ip is not a strict impact-parameter rep- 
resentation of the scattering amplitude. This is because 
of the dependence of S'(b) on additional kinematic vari- 
ables like K or q which we do not show explicitly and 
the angle-dependent factor K — k cos(6'/2) in front of the 
impact-parameter integral: in a genuine impact-parameter 
representation all dependence on the scattering angle 9 
should only reside in the factor exp(— iq • b) 27J. 

The main features of these representations are func- 
tional integration over velocities without boundary con- 
ditions and the use of "phantom" degrees of freedom to 
get rid of explicit phases which would diverge in the limit 
of large scattering times. Two versions exist which are 
distinguished by the reference path along which the par- 
ticle dominantly travels and the dimensionality d of the 
"anti- velocity" w(i) which is needed to achieve the can- 
cellation : 

{+00 
I j d<y [v2(t)-w2(t)] 
— oo 

X exp| i x(b, V, w] | (2.4a) 

+00 

X(b,v,w] = - j dtV{x,,f{t) + 

Equant 

(i, V, w] ) .(2.4b) 

—00 

In the first case the reference path is a straight-line path 
along the mean momentum 



K 



'(t) = h+-t 
m 



and the quantum fluctuations are given by 



(d=3) 
quant 



{t,v,w] = x„(t) - x„(0) 



(2.5) 



(2.6) 



where 
x„(t) 



+ 00 



dt'sgn{t-t')w{t'), i,„(i) = v(t) (2.7) 



and sgn(a;) — 20{x) — 1 is the sign-function. We will call 
that the "aikonal" representation because it gives rise to 
the eikonal approximation of Abarbanel & Itzykson (AI) 
[28] if the quantum fluctuations are neglected altogether. 



^ The path integrals are normalized such that S{h) = 1 for 
zero potential. Our notation indicates that x is a function of 
b but a functional of v(f) and w(t). Similarly for Xquant. 
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In the second case the anti- velocity is only 1-dimensional 
and the reference path is a ray along the initial momentum 
for i < and along the final momentum for t > Q 



^lti'\t) = b + [k. ei-t) + kf e{t) 

m Zm 

and the quantum fluctuations are given by 

Xqulilt V, w] = x„(t) - x^„(0) - K a; 11^,(0) 



(2., 



(2.9) 



We will call that the "ray" representation in the following. 

For further details we refer to ref. [T7]- Here we just 
note that systematic eikonal-like expansions can be ob- 
tained by Taylor-expanding the potential around the ref- 
erence path 



X(b,v,v^f] = - 



1 



+ 00 



dt 



F(Xref) + Xquant ' VV^(Xrcf ) 



~r Equant )i (Xquant)j didjV {-X-rd) + . . . (2.10) 

and performing successively the functional (Gaussian) in- 
tegrations over velocity and anti- velocity. This is because 
a simple scaling argument shows that each quantum fluc- 
tuation Xquant IS supprcsscd by a power of 1/Vk in the 
case of the "aikonal" representation or 1/Vk in the "ray" 
representation. Therefore at high energy/small scattering 
angle the geometrical (classical) picture of scattering is 
dominant. 



3 Variational calculation 

3.1 The Feynman-Jensen variational principle 

Being highly nonlinear in velocity and anti- velocity vari- 
ables the path integrals in eq. (|2.4ap cannot be performed 
analytically in general 0. Numerical methods or approxi- 
mations are then necessary. For weak potentials, for ex- 
ample, one may expand in powers of the potential and it 
has been shown that the Born series for the T-matrix is 
obtained in all orders 17 . At high energy and small scat- 
tering angles eikonal approximations are useful. Because a 
variational approach neither requires weak interaction nor 
high energy, forward scattering it is widely used in atomic 
and molecular physics, mostly in the form of Kohn's vari- 
ational principle |7j. Schwinger's functional is also station- 
ary against variation of trial wave functions [B] but more 
difficult to use in practice. 



It would be interesting to derive the few exact expres- 
sions for scattering amplitudes of a local potential which are 
available in the quantum-mechanical literature, e.g. for the 
Coulomb potential. 



The path-integral representation (I2.4ap immediately 
suggests another variational approximation for the action 



A[v, w] 



dt 



V[ Xrof(t) + Xquant (i,V,w] 



(3.1) 



viz. the Feynman-Jensen variational principle. For a pos- 
itive weight function this principle employs the convexity 
of exp(— x) to obtain the inequality 



(3.2) 



but in real time (in which scattering occurs) one only has 
stationarity: 



5(b) 



where 



V^vV^we 



/ V^vV^w exp {iAt +iA- lAt) 
/ V^vV^w exp {iAt) 

St{h) ■ e'<^-^>t , AA = A-At (3.3) 



_ jVWw AA exp(zA) , . 
\^^)t- Jps^pdy; exp(zA) ^ ' 

is the average of the difference between the full action 
A and the trial action At weighted with the oscillating 
factor exp(iylt). Note that both Kohn's and Schwinger's 
variational principle also are only stationary when applied 
to the full T-matrix. This seems inevitable when trying to 
estimate a complex quantity by variational means; only 
for real quantities, like a scattering length, a minimum 
principle is available (see Chapt. 11.3.4 in ref. @]). 

3.2 Variational ansatz and equations 

As usual in variational calculations the outcome crucially 
depends on the test functions/actions which must lead to 
expressions which may be evaluated safely. Unfortunately 
in the path-integral formalism one is restricted to trial ac- 
tions which are at most quadratic in the dynamical vari- 
ables so that the various path integrals and averages can 
be worked out analytically. This restricts somewhat the 
utility of this approach which - on the other hand - is so 
general that it can be applied not only for the scattering 
of a single particle in quantum mechanics but also in the 
many-body case or even in field theory. In the present case 
it means that our trial action At may contain linear and at 
most quadratic terms in the velocity and the anti- velocity. 
So we may take 



A[v, w] 



+ 00 



dtdt' — v{t)Ayy{t,t')v{t') 

+ w{t) A„„(t,i')w(t') 
+ yr{t) A^^{t,t')w{t') 

-f w(i) A^^{t,t')w{t') 



-\-oo 



dt 



B„(t)-v(^)+B„(t)-w(t) (3.5) 
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since constant terms cancel out in the Feynman- Jensen 
variational principle. Special cases are the free action and 
the ansatz studied in ref. [25] where only the B-terms were 
varied while the quadratic part was left as in the free case. 

At this stage it is useful to have a look at the high- 
energy expansion in eq. (|2.10p : It tells us that the func- 
tions Ayy , Auju; , Aym , Ayj y havc to bc anisotropic in order 
to describe the high-energy expansion up to and including 
terms of order (in the "aikonal" representation) or 

in the "ray" representation. Thus the various varia- 
tional functions Ayy (i, t') . . . actually are 3x3 matrices in 
cartesian space if the most general quadratic trial action is 
considered. Similar anisotropic trial actions have already 
been considered in a variational description of world-line 
scalar field theory [5J and shown to give considerable im- 
provement. A slight complication here is the presence of 
the anti-velocity which, however, can be elegantly han- 
dled by grouping it together with the velocity to form a 
(3 -|- d)-dimcnsional "super-vector" 



V(3-3) ^ / V 



or 



V 

W|| 



(3.6) 



SO that the trial action can be written succinctly as 

771/ 

At = — V-AV + B-V . 
2 



(3.7) 



Here IB is a (3 -I- (i)-dimensional vector made up of the 
variational functions B(t) 



B„ 



(3.8) 



and A a (3 -|- d) x (3 -t- (i)-dimensional symmetric matrix 
formed by the variational 3x3 matrix functions A{t, t') 



(3.9) 



In eq. p.7p we employ a nomenclature where also the inte- 
gration over continuous times is treated like a summation 
over identical indices. See app. El for a more detailed ac- 
count of our conventions. This allows to evaluate the var- 
ious path integrals and averages in an efficient way and to 
derive the variational equations easily as detailed in app. 

m 

Here we just collect the final results: after variation the 
stationary value of the impact-parameter iS-matrix reads 



5var(b) = exp [t{Xo+Xi)-f2] (3.10) 



where 



-\-oc 



(3.11) 



J —oc 

1 r+°° 

Xi^- dtdt' d,v,{t){p{t)) {So),^ {t,t') 

xd,V^^t,){p{t')) = ^jJ SoJa (3.12) 



are "phases" and 
1, 



(3.13) 



is the "fluctuation term" (or an imaginary phase) which 
arises from the functional determinant due to the quadratic 
term in our ansatz. "Tr" denotes the trace both in continu- 
ous and discrete variables and "T" the transpose of a vec- 
tor or a matrix. Actually, the terminology used above is a 
little bit misleading as it turns out that all these quantities 
become complex. This is because they are not determined 
by the real potential and its derivatives but by its Gaus- 
sian transforms with a complex width a{t) (see below). 
This entails that different branches of the multi-valued 
logarithmic function "Ln" may be needed depending on 
the appropriate sign of the square root of the complex 
functional determinant. These subtleties of the complex 
Gaussian integral which require a "branch tracking" are 
discussed in app. ID. 41 

The Gaussian transform of the potential is most con- 
veniently defined in momentum space as 



KT(t)(p) Vip) exp 
= V{p) exp 



1 



p^(7(i)p 



;Pia,j{t)pj 



(3.14) 



and we can form the Jacobian (the vector of derivatives) 

(J,), d, V, (3.15) 

and the Hessian (the matrix of second derivatives) 

[R,)^^ := d,d, K . (3.16) 

They are functions of a trajectory p^^^(t) (we will omit 
the subscript "var" in the following) which is determined 
by the following equation of motion 

/ + 00 
dt'Z'o(t,t')Ja(t')(Pvar(^')) ■ 
-oo 

(3.17) 

Here we encounter the (3 x 3)-matrix Uq which in the 
"aikonal" representation takes the form 



4^-3) (t,i') 



1 

2m 



1 

\t~t'\ I 10 
1 



(3.18) 



whereas in the "ray" representation it reads 



, / d{t,t') 
4'''\t,t') ^ -T^ \ d{t,t') I (3.19a) 
2™ \ \t~t'\ 

with 

d{t,t') = \t-t'\ - \t\ - \t'\ (3.19b) 

if the mean momentum K is chosen to be along the 3- 
direction. Finally, it is found that the time-dependent width 
of the Gaussian transform 



a{t) = iS{t,t) 



(3.20) 
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is proportional to the diagonal part of a matrix S which 
fulfills a Lippmann-Schwinger-like equation 



U — Uq + Sn Her S . 



(3.21) 



In contrast to the familiar Lippmann-Schwinger equation 
in scattering theory this is a nonlinear equation since the 
Gaussian width of the potential and its Hessian depends 
on the diagonal part of S. Note that eq. p.2ip holds for 
both representations and that all dependence on the phan- 
tom degrees of freedom has now disappeared. In particu- 
lar, no divergences are encountered for large scattering 
times - the job for which the anti- velocity was introduced 
in the beginning. We find it also remarkable that in the 
final outcome the variational functions A, B whose physi- 
cal interpretation is not immediately evident are replaced 
by the trajectory p and the Green function E which have 
a clear classical and quantum-mechanical meaning. 



Eqs. (PTfl) . ((X^ and ([5^ form a system of coupled 
variational equations which have to be solved. Once that 
is done we may insert the solutions into the variational 
functional to get the impact-parameter iS-matrix (I3.10p . 



3.3 Correction by the second cumulant 

It is possible to calculate systematic corrections to the 
variational result (or improving the variational principle) 
by realizing that the Feynman- Jensen variational principle 
is the first term of a cumulant expansion 



exp 



iAi + ^ ^2 + 



where 



Ar = {AA)^ 



(3.23) 



is just what enters the Feynman- Jensen variational prin- 
ciple and 



(3.24) 



the first correction to it. So the variational functional be- 
comes 



vai+cum 



(b) 



exp 



(3.25) 



In principle one could vary the full functional with respect 
to the variational functions/parameters but we follow the 
standard practice to add the second cumulant as a cor- 
rection |25) with the parameters fixed by the Feynman- 
Jensen variational principle, i. e. the first cumulant. This 
allows to use the variational equations for these parame- 
ters and leads to a considerable simplification of the final 
expression. 

App. Ogives the result of evaluating the second cumu- 
lant in our case: 

A2 = {x'),-{xft~^Jl ^J- + lT^H^H,f (3-26) 



where 



dti dt2 



X exp-^ 



,T „T^ ( S{ti,ti) S{h,t2)\ ( Pi 



(pr,pD 



+ i 




( T 
(Pl .P 



(3.27) 



involves a double Gaussian transform of the squared po- 
tential. All other quantities have been defined and cal- 
culated before in both representations. It is easy to see 
that AA is of first (and higher) order in the potential. 
By construction the n*'^ cumulant then contains contribu- 
tions of ©(y") and higher. This implies that the second 
Born approximation to the T-matrix is fully included in 
our calculation when the second cumulant is added. 



3.4 A special case: The linear ansatz 

A less general variational ansatz has been made in ref. 
[22j by only allowing the linear terms in the trial action 
to vary whereas the quadratic part was fixed to be the 
free action. In the present nomenclature this amounts to 
setting A = (T3. Inspecting the variational solution for A 
we see that the results of ref. [H] should be recovered 
by setting the Hessian of the potential to zero: TJo- — >■ 0. 
Indeed, then the solution of the Lippmann-Schwinger eq. 
p. 211) simply is S ~ Sq and the Gaussian width becomes 



(3.22) a^jit) 



it 





\ (1 - 5.3) 



"aikonal" /o ooa 
"ray" . ^^■^^> 



Thus in the "aikonal" case (with a 3-dimensional anti- 
velocity) there is no Gaussian transform of the potential 
and the variational equation for the trajectory simply is 



Pit) 



b + — i 



1 

2m 



+ 00 



dt' \t-t'\^V{p{t')) (3.29) 



which, after differentiating twice, is just Newton's law for 
the classical motion in the potential: 



(3.30) 



However, the boundary conditions encoded in the integral 
eq. p. 291) are unusual: using \t — t'\ — > \t\ — sgn{t)t' for 
large |t| one finds from eq. p.29p 



P{t) 



i— !-±oc 



K 

— T 



-|-'DO 



dt 



2m 



+00 



t + h± 



dt't''^'^ 



2m 



terms which vanish for |t| — )■ 00 ,(3.31) 



so that 



hm {piT)+p{-T)-T[p{T)-p{-T)]} = 2h (3.32a) 

lim {p(r) + p(-T)} = 2— (3.32b) 

T->oc m 
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It is the great advantage of the integration over velocities 
that the boundary conditions are exactly built in and we 
do not have to worry about them during (approximate) 
evaluation of the path integral. Note that the boundary 
conditions (I3.32ap and (|3.32bp are just those needed to 
convert an integral equation with the kernel \t — t'\ into 
a differential equation 29j. 

For the linear ansatz in the "aikonal" representation 
the trajectories are real and consequently the phases Xq 
and Xi are also real; an imaginary phase, i.e. a real part in 
InS'(b) only develops due the second cumulant. It is easily 
seen that the general expression p.27p for A2 reduces to 

the eq. (4.84) in ref. gg, if T is replac ed by 17^^"^^ As 
usual one finds from Newton's equation p.30p the conser- 
vation law mjy^{t)/2 -f V{p{t)) — E — const. Inserting 
the asymptotic behaviour p. 311) and using the symmetry 
of the solutions under t ^ —t which forces J^^dt' J(t') 
to be perpendicular to the mean momentum K, the con- 
served quantity is 



K2 1 

E = — + — 
2m 8m 



4-00 



di' J(t') 



(3.33) 



which means that this "energy" depends on the impact 
parameter b. This is because 5'(b) — S'(b, K) does not 
contain any information on the actual scattering energy 

2m 2m 8m 

However the variational principle is "clever" enough to 
mimic the missing term as good as possible: 



(3.34) 



di' Vbl/(p(t')) 



(3.35) 



is what one would obtain if the impact-parameter integral 
over exp(iAro) is evaluated by a saddle-point approxima- 
tion. 

Things are slightly different in the "ray" representa- 
tion: First, the potential is replaced by its Gaussian trans- 
form with an anisotropic width given in eq. p.28p which 
renders it complex from the very beginning. Second, the 
anisotropic kernel p.l9ap leads to an equation of motion 
where the particle experiences a kick at t = and the 
"energy" of that motion is not conserved anymore during 
the scattering process but only asymptotically (see fig. 1 
in ref. [22]). Complex trajectories are frequently encoun- 
tered in semi-classical approximations of quantum motion 
in regions of space which are forbidden classically. In our 
variational approximations, however, the main reason for 
this behaviour (which shows up even above any potential 
hill) is unitarity: ^(b) cannot be unimodular |18) . 

We do not dwell on several other interesting properties 
of our variational approximation with a linear ansatz but 
refer to ref. [22l for more details. 



3.5 High-energy expansion 



energy expansion (I2.10p . To show how that comes out we 
will take here for simplicity the "aikonal" representation 
as the algebra in the "ray" representation is more involved. 

First, we switch to distances instead of times by defin- 
ing 

K m , ^ 

z —t ^ t = —z. (3.36) 



m 



K 



We then see that 



.(3-3) 



(z,z') 



2K 



(3.37) 



is suppressed for large energies and forward scattering an- 
gles (recall K = fccos(6'/2)) since the arguments z,z' are 
bounded by the range of the potential. We thus may ex- 
pand in powers of Sq and obtain from the Lippmann- 
Schwinger-like eq. p.2ip 

S = Sq + SoH^So + ■■■ ■ (3.38) 
Eq. p.20p then gives for the Gaussian width 



CTy(z) 



dz'\z-z'\' ^{H,^^,^)^^{z') + . 



my 

(3.39) 

as the first-order term vanishes in the "aikonal" represen- 
tation (this is not the case in the "ray" representation !). 
Therefore the Gaussian width becomes also small under 
these kinematical conditions (note that m/ K doesn't need 
to be small) so that may be replaced by H up to order 
From the definition p.l4p of the Gaussian trans- 
form wc find 

V;(.)(p(z)) = V{p{z)) + ^o,,{z)HM^.)) + ... (3.40) 



and iterating the equation of motion p.l7p one obtains 

p = x,ef + + HSq J + . . . . (3.41) 

Here and in the following the potential, the Jacobian and 
the Hessian are all evaluated with the straight-line refer- 
ence path Xicf (z) = b -|- Kz . Inserting these expressions 
into eqs. p. lip - p.l3p . and expanding up to order 
we get in our short-hand notation 

Xo ^ XaI - J^^oJ - p^SoHSoJ - ^Tr (Z'oiI)'(3.42) 



Xi 
Q 

where 



(0) 

Xai 



1 



-J SqJ + J SqHSqJ 



(3.43) 



+ 00 



dtV{y.,^t{t)) = -— I dzV[h + Kz 



+ 00 



(3.45) 

is the standard eikonal phase from Abarbanel & Itzykson 
[2H1. We see that 



By construction the variational solution of the quadratic 
-|- linear variational ansatz p.7p must contain the high- 



X(i + Xi 



(2)^. (2) 



O 
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where in our condensed notation 



3.6 A Feynman-Hellman theorem 



Xai 

,(2) 



1 ^ 

1 



(2) 

Xai 



(3.48) 



are exactly the first- and second-order corrections to the 
leading "aikonal" result as given in eqs. (4.97), (4.103) 
and (4.110) of ref. [H]. Note that the contribution (P^IH) 
from the fluctuation term brings the imaginary part into 
full agreement with what is needed for unitarity in the 
high-energy limit: 



(2) 
AI 



,(2) 
^AI 



o 



1 

7^3 



,(2) 
^AI 



o 



1 

' (3.49) 

How is that high-energy behaviour changed by includ- 
ing the second cumulant? To answer this question we con- 
sider eq. p. 271) and decompose 



1 / <j{ti) 
I \ a{t2 

S{tiM) 



(3.50) 



into a diagonal and a nondiagonal part. The diagonal part 
is seen to convert the potentials into Gaussian transformed 
potentials while the nondiagonal part gives new contribu- 
tions when the exponential is expanded: 



+ 00 



dtidt2 {d2iS{ti,t2)ijdij) 

xv,(^t,){p{h))v^(t,){p{t2)) + 



.(3.51) 



It is seen that the the first 3 terms on the r.h.s. are com- 
pletely cancelled when forming the second cumulant in eq. 
p.26p and that in the high-energy limit the leading con- 
tribution comes from the last term involving 3 powers of 
the matrix S. In that limit we may safely replace it by the 
zeroth-order matrix Sq so that the second cumulant gives 
an imaginary (phase) correction at 0{1/K^). It is inter- 
esting that in the linear approximation the last term in 
eq. p.26p is absent and that the contribution — A2/2 from 

the second cumulant then correctly supplies the real term 
(2) 

uiaj- In contrast, the more general, anisotropic quadratic 
-|- linear trial action already gives variational results cor- 
rect up to and including 0{1/K'^) as expected from the 
high-energy expansion. Similar results hold in the "ray" 
representation except that here the expansion parame- 
ter is the wave number k and not the mean momentum 
K = kcos{e/2). 



* To show that one has to convert to an impact-parameter 
representation of the scattering matrix for which the K = 
fc cos(0/2)-factors in the "aikonal" representation give addi- 
tional contributions. As they are the same as worked out in 
ref. [22] we do not have to consider them here. 



Some further properties of the variational approximation 
for the impact-parameter 5-matrix (13.101) can be obtained 
by employing the variational equations in a particular way. 
This has already been used in a variational approximation 
to the relativistic bound-state problem |30] and is just a 
variant of the well-known Feynman-Hellmann theorem in 
quantum mechanics (see, e.g. refs. [31,32 ). 

Suppose there is some parameter A in the potential 
or in the reference path. As the variational functional 
also depends on it implicitly via the variational param- 
eters/functions Avar and Bvar ouc may perform the differ- 
entiation with respect to this parameter by means of the 
chain rule 



91nSV 
dX 



d In Svai 






dX 


r+00 


, d 






"V 

^ — 00 


''dX 


p-\-oo 


, d 






J —00 


''dX 



d In gyar 

(tMt)) 
(Pit)) 



gvar 



dX 



a(t) 



(3.52) 



since the first two terms vanish identically for the vari- 
ational solutions Avar and Bvar- This leads to a simple 
result although these solutions may be very complicated 
functions of the parameter A. 

As an example take the dependence of the impact- 
parameter iS-matrix on the mean momentum K in the 
"aikonal" representation. The sole explicit dependence of 
the variational functional resides in the reference path 
Xrcf(i) = b + Kt/m. Therefore we immediately have the 
relation 



dK 



m 



dttVV,(t){p{t)). (3.53) 



Similarly one obtains 

d In S'var 



db 



ih 



+ 00 



dtVV,(,){p{t)) (3.54) 



which also holds for the "ray" representation. Finally the 
dependence on the strength Vq of the potential is given by 
the simple expression 



ainSy 



dtV,(^t){p{t)) = iXo. (3.55) 



These relations may be used to derive or check high-energy 
or weak-coupling expansions by expanding the r.h.s. in 
powers of the particular parameter and integrating term 
by term. 



^ Here it is important that only Xq is assumed to depend 
explicitly on A whereas Xi,f} do not. See eqs. ()B.2|l . (|B.3|) . 
(|B.4|) in app. [Bl We use the short-hand notation of app. [X] 
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For example, using the first two terms of the high- 
energy expansion (j3.41l) we have in the "aikonal" repre- 
sentation 



d 



1 



K.VK.(,)(p(z)) = —yv{^,,i{z))- — yv{^,,i{z)) 

dz' \Z - z'\ VF(Xref (z')) + 0{1/K'^) 1 (3.56) 



since the Gaussian width does not contribute in this or- 
der. Substituting this into the Feynman-Hellmann relation 
(j3.53p (expressed in the z-variable to make Xjof (z) inde- 
pendent of A') and integrating with respect to K gives the 

correct resuh In^iar"^^ = ix^^] + ix^^] + 0{l/ K^) if an 
integration by parts in z is performed. 

In aPD. ID.61 we use the dependence (j3.55D on the cou- 
pling strength as a nontrivial check for our numerical pro- 
cedures to solve the variational equations. 



3.7 Semi-classical expansion 

It is of some interest to distinguish the present variational 
approach from the usual semi-classical approximation of 
the path-integral propagator obtained by the stationary 
phase approximation (see, e.g. ref. [33], Chapt. 13 or ref. 
[33], Chapt. 4). To do that we have to restore the appro- 
priate factors in the path-integral representation (|2.4ap . 
This is easily done by remembering that the weight in the 
path integral is always the exponential of i x the action 
divided by h and that the velocities in the reference path 
are momenta {— hx wavenumbers) divided by the mass 
of the particle. To account for that we simply have to 
substitute 



V 



1 



(3.57) 
(3.58) 



everywhere. From the explicit expressions (I3.18|) and p.l9a[l 
we then see that 

= 0{h) (3.59) 



and from the Lippmann-Schwinger-like equation p.2ip 
that 

S = 0{h) + higher orders in h (3.60) 



since eq. p.58p also implies the multiplication of the Ja- 
cobian and Hessian by 1/h. In other words: the Gaus- 
sian width (7{t) = S{t,t) is a pure quantum effect which 
vanishes in the classical approximation. The "higher or- 
ders in h " in eq. p.60p are generated by the nonlinearity 
of the Lippmann-Schwinger-like equation via the Gaus- 
sian transform of the Hessian . Similarly the Gaussian 
transform of the Jacobian J„ in eq. p.l7p gives rise to 
/i-dependent terms in the trajectory 



Pvar(0 - Pcl.ss{t)+0{h) + 



(3.61) 



where Pciass(^) the trajectory for zero Gaussian width. 
By the same arguments we see from eqs. (j3.1ip . (|3.12D 



that the phases 



o 



whereas from eq. p.l3p one finds that 
n = 0{h°) + ... . 



(3.62) 



(3.63) 



Finally, applying the scaling laws (|3.58p and p.60p to eqs. 
dSmj), ([3^ we deduce that 



A, 



0{h) 



(3.64) 



is a genuine quantum correction. We thus see that our 
variational approach is not equivalent to a semi-classical 
approximation as it contains (an infinite number of) h- 
dependent terms. This is mostly due to the quadratic term 
in our general variational ansatz p.7p which gives rise to 
an interaction-dependent Gaussian width: the more re- 
stricted linear ansatz studied in ref. |22| leads to a van- 
ishing Gaussian width in the "aikonal" (3-3) representa- 
tion; only in the "ray" (3-1) representation one obtains a 
nonzero width which generates some quantum corrections. 



4 Numerical results 

We compare our results with those of a (basically) exact 
partial- wave calculation of scattering from a Gaussian po- 
tential 

V{r) ^ Vbe-""^' = Voe""-'/^' (4.1) 



with parameters 



2mVoR^ 



kR^4 



(4.2) 



This was considered as a test case for the systematic eikonal 
expansion by Wallace in ref. [TH] and (with Vq = —41.6 
MeV and R = 1 fm) used to describe phenomenologically 
the scattering of a-particles from a-particles at an energy 
of 166 MeV. We also choose it because for this potential 
a persistent failure of the eikonal expansion to describe 
the scattering at larger angles was observed: the scatter- 
ing angle where the eikonal amplitude started to deviate 
appreciably increased only slightly when higher-order cor- 
rections were included (see fig. 6 in ref. |18j). 

Thus this particular potential is a good litmus test for 
checking any approximate description of high-energy scat- 
tering. However, we do not use 



Aa 

a 



da /df2\g 



(i(T/(ii7|oxact 



da /dn\c 



(4.3) 



as the usual measure to gauge agreement/disagreement 
with the exact result but 



(4.4) 







/approx /exact 


/ 




/exact 
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Here 

f{o) = -^r.^f (4.5) 

is the scattering amplitude and 

% = 1/(^)1^ (4-6) 



Table 1. Exact partial- wave amplitude and differential cross 
section for scattering from the Gaussian potential (|4.ip with 
2mVoR^ ~ —4 at kR — 4. Amplitude and cross section are 
given in units of R and R^, respectively, where R is the range 
of the potential. The number in parenthesis indicates the power 
of ten by which the numerical value has to be multiplied, e.g. 
5.6472 (-6) = 5.6472 • 10"^ 



the differential cross section. By construction the quan- 
tity in eq. (j4.4p also measures the phase deviation of the 
scattering amplitude /approx from the exact amplitude. In- 
deed, if 

/exact = re'-^ and /app,ox = (r + /Ir) e^(^+^^) (4.7) 
then one finds to first order in the deviations Ar, A<f> 



er] {Ref)/R (Im/)/i? 



f 



^{A^f + iAr/rf (4.8) 



and 



Act _ 2^*" 
a r 



which is independent of the phase error A(j). As 



Aa 


< 2 


^/ 


a 




/ 



(4.9) 



(4.10) 



one may even encounter cases where Aa/a = although 
7^ because the approximate scattering amplitude 
has the correct modulus but the wrong phase. So the com- 
mon practice of comparing just differential cross sections 
can be misleading. 

Table [1] gives the values of the exact calculation ob- 
tained by integrating the Schrodinger equation for each 
partial wave up to a radius where the solution could be 
matched to a linear combination of free spherical waves. 
By varying this radius as well as the step size for the 
numerical integration and the number of partial waves 
retained we made sure that the numerical values given 
in this table are at least accurate to the number of dig- 
its given. We display scattering amplitude and differential 
cross section in units of powers of R since the range of the 
potential sets the length scale. 



4.1 Variational results 

We have evaluated the coupled variational equations by 
an iterative scheme on a grid in time (more precisely: in 
distances which the particle travels while in the range of 
the potential) and then evaluated the phases Xo,Xi, the 
fiuctuation term f2 and the second cumulant A2. Some 
numerical details are given in app. |D] 

In table[2]we list the relative deviation of the scattering 
amplitude for the various approximations from the exact 
partial- wave result as a function of the scattering angle 9. 
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For the different variational approximations we use the 
following naming scheme: 

B 3d : variational ansatz with linear term B 
only ED , 
Be 3d : variational ansatz with linear term B 

-|- second cumulant [H] , 
AB 3d : variational ansatz with quadratic 

-I- linear terms A and B , 
ABc 3d : variational ansatz with quadr. -I- linear 
terms A and B -(- second cumulant , 

(4.11) 

where d = 3 ("aikonal") or d = 1 ("ray"). Note that we 
list 



10^ 



/ 



(4.12) 



which is necessary to show the rather small deviations of 
the scattering amplitude - any discrepancies in the cross 
sections wouldn't be seen on a logarithmic scale except at 
large scattering angles. 

Two independent programs using different subroutines 
and (partially) different techniques have been developed 
and run to arrive at these results. Many checks have been 
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Table 2. The relative deviation (in units of 10~'') of the scattering amplitude from the exact value for scattering from the 
Gaussian potential (|4.1[) with 2mVoR^ = —4 at kR = 4. The variational approximations are labeled according to the scheme 
()4.11|1 . The numbers in parenthesis give the estimated numerical error of the calculations in units of the last digit. Entries are 
stopped if the deviation exceeds 10 %. 



e [°] B 33 B 31 Be 33 Be 31 AB 33 AB 31 ABc 33 ABc 31 
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performed (see app. D) including a test of the Feynman- 
Hellmann theorem for the interaction strength which is a 
highly nontrivial test of how well the variational equations 
are fulfilled. In addition, the various accuracy parameters 
have been varied to ensure stability of the numerical out- 
come. For the approximations including the second cumu- 
lant which demand particular care we also include in table 
[2] an approximate numerical error estimated from the dif- 
ferent results from both programs and from variation of 
the integration parameters. 



A systematic improvement is seen when quadratic terms 
and/or the correction by the second cumulant are included. 
This is shown in fig. [T] were the deviations of the different 
results for the scattering amplitude are plotted on a linear 
scale. We also have included the best result from the varia- 
tional approximation with a linear ansatz plus the second 
cumulant (Be 31). It is seen that including the second cu- 
mulant gives much improved results so that the variational 
approximations labelled Be 31 and ABc 31 give the best 
overall results. However, at present we do not fully under- 
stand why at larger scattering angles the (much simpler) 
scheme Be 31 outperforms the more involved scheme ABc 
31 which includes the quantum- mechanical spreading ef- 
fects. As we simply have added the second cumulant on 
top of the variational results and not used it in the varia- 
tional optimization there is, of course, no guarantee that 
ABc 31 does always better than Be 31 whereas AB 3d has 



always to be better or at least as good as B 3dE|. It should 
also be noted that the evaluation of the second cumulant 
at large scattering angles where the amplitude is down by 
many orders of magnitude is a very demanding numeri- 
cal task and the numerical errors given in table [5] may be 
underestimated . 

Nevertheless, the full "ray" approximation plus the 
second cumulant (ABc 31), for example, provides a very 
good approximation even at larger scattering angles: at 
6 — 90° the complex scattering amplitude is reproduced 
within 9 % and the cross section to better than 2 % al- 
though the latter one is down by 5 orders of magnitude 
from the value taken in forward direction. 

Only at large scattering angles the deviations from the 
exact partial-wave result again start to grow and become 
so large that they even can be seen on the usual logarith- 
mic plot for the cross section in fig. [H 



4.2 Comparison with other work 

Combining the eikonal approximation with the second- 
order Born approximation Chen has proposed several ap- 
proximations which should also work at larger scattering 
angles and which we denote by "Chen 1" [3S], "Chen 2" 
[SS] and "Chen 3" [37. 



Otherwise the variational principle would have chosen the 
free A as the best solution. 
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6 (degrees) 

Fig. 1. (Color online) The relative deviation \Af / f\ of the variational amplitudes from the exact scattering amplitude as 
function of the scattering angle 9. The naming scheme (|4.11|) is used for the various approximations. 



The relative deviation of his approximate amplitudes 
from the exact value has been calculated by direct nu- 
merical integration and is displayed in table [3] and fig. [31 
It is seen that these approximations are clearly inferior 
to the present variational approximations although they 
use some analytic properties specific to the potential un- 
der consideration. In contrast, the variational approach 
requires as input just the (Gaussian transform of the) po- 
tential and thus can be applied more generally. 

We also include the results for Wallace's systematic 
eikonal expansion Q including real and imaginary phases 
up to second and third order in 1/k. This makes for a rel- 
evant comparison with our variational results which have 
been shown to be correct up to second order in the high- 
energy expansion but contain many other higher-order 
terms. 

Table |3] and fig. [3] show that these systematic expan- 
sions are very good in forward direction but fail at higher 
scattering angles. Again the variational approximations 
are clearly better in this kinematical regime. 



For a Gaussian potential the explicit expressions are given 
in table I of ref. [H]. 



5 Summary and outlook 

We have applied the Feynman- Jensen variational princi- 
ple to two variants of a new path-integral representation of 
the T-matrix in potential scattering using the most gen- 
eral trial action which contains quadratic and linear terms 
both in velocity and anti-velocity. Merging these variables 
into a "super-velocity" the variational functional was eas- 
ily derived and the variational equations immediately fol- 
lowed by demanding stationarity of that functional. As ex- 
pected the anti- velocity ("phantom") degrees of freedom 
prevented any divergence to appear when the scattering 
time was sent to infinity as it should be when the initial 
and final states are detected in the scattering process. 

We have found that this quadratic -I- linear ansatz de- 
scribes the scattering in terms of a (nearly) classical tra- 
jectory p{t) with peculiar boundary conditions and a ma- 
trix operator S which is the solution of a one-dimensional 
Lippmann-Schwinger-like equation. This reflects the clas- 
sical as well as the quantum- mechanical ("wave spread- 
ing" ) aspects of the process dominant in the high-energy, 
forward scattering and low-energy scattering regime, re- 
spectively. Unlike the previous linear ansatz |22| the present 
variational approach immediately leads to complex phases 
in the impact-parameter iS-matrix as required by unitar- 
ity and a fluctuation term coming from the corresponding 
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Fig. 2. (Color online) The exact cross section for scattering from a Gaussian potential with strength 2mVoR^ = —4 at kR — 4 
and the best variational cross sections including the second cumulant. 



functional determinant. We were also able to evaluate the 
first correction (the second cumulant) in a systematic cu- 
mulant expansion. 

As we solved the nonlinear, coupled variational equa- 
tions iteratively (both in analytic and numerical evalu- 
ations) the present work is mostly devoted to the high- 
energy case with "wave-spreading" corrections becoming 
more important at backward scattering angles. Indeed we 
have seen better and better agreement with exact partial- 
wave calculations of scattering from a Gaussian potential 
as we step from the linear trial action to the quadratic -I- 
linear and from the version which uses a straight-line ref- 
erence path and a 3-dimensional anti- velocity ("aikonal" 
or "3-3" ) to the version with a ray-like reference path and 
a 1-dimensional anti- velocity ("ray" or "3-1"). An even 
more dramatic improvement is obtained if the correction 
by the second cumulant is added because then the second- 
order Born approximation is contained in the variational 
result - an observation also made by Chen (refs. [3S] - [37] 
) . Instead of plotting the approximate cross sections on a 
logarithmic scale as is usually done, we give the numerical 
values for the relative deviation of the approximate scat- 
tering amplitudes from the exact partial- wave result. This 
is done for a particular benchmark case: scattering from a 
Gaussian potential with fixed parameters (|4.2p and allows 
a quantitative comparison of the different approaches for 



a case where large-angle scattering has been notoriously 
difficult to describe. 

We also have shown analytically that in a high-energy 
expansion terms up to inverse second order in k or K = 
kcos{9/2) for "ray" or "aikonal" representation are cor- 
rectly reproduced. Thus our variational approximation -|- 
cumulant correction is correct up to next-to-next-to-leading 
order in a systematic eikonal expansion. 

However, there is no reason to limit the variational 
approach to high energies only since we expect it to be 
useful at lower energies also. Indeed, it is just one of the 
inherent advantages of a variational description over an 
expansion in some small parameter that the variational 
principle stretches the parameters/functions in the trial 
ansatz in such a way as to make it "optimal" even under 
unfavorable kinematic conditions. Unfortunately, we were 
not able to corroborate this expectation either by analyt- 
ical or numerical methods. Obviously, a better strategy 
to solve the highly nonlinear variational equations is re- 
quired. 

This is one of the tasks left for the future. The other 
one is, of course, to apply this approximation scheme to 
more realistic situations, e.g. many-body scattering. Since 
the case of a spherically symmetric potential V{r) was 
only chosen to compare with a (numerically) exactly solv- 
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Table 3. Relative deviation (in units of 10^*) of the scattering amplitude from the exact value for various approximations by 
Chen [3511361157] and Wallace [TS]. In the latter case the subscripts denote the order of the eikonal expansion which is included. 
The potential is again a Gaussian one with the same parameters as in table [J] 



e [°] Chen 1 Chen 2 Chen 3 Wallace2 Wallaces 






217 


5 


40 


4 


< 1 


5 


220 


6 


40 


4 


1 


10 


231 


9 


42 


5 


1 


15 


250 


15 


45 


7 


1 


20 


279 


25 


49 


9 


1 


25 


320 


42 


55 


9 


2 


30 


378 


67 


63 


8 


4 


35 


456 


108 


72 


3 


9 


40 


560 


170 


84 


17 


15 


45 


693 


266 


98 


53 


24 


50 


844 


406 


115 


122 


35 


55 


971 


585 


142 


232 


47 


60 


1001 


767 


194 


376 


66 


65 


891 


907 


281 


528 


110 


70 


672 


999 


405 


678 


192 


75 


400 


1078 


563 


844 


317 


80 


244 


1190 


765 


1065 


499 


85 


555 


1390 


1024 


1392 


763 


90 


1082 


1739 


1357 


1898 


1150 



able model, there is no difficulty to apply our variational 
scheme to general potentials and scattering from few-body 
targets. A good step into that direction would be to com- 
pare with the Faddeev results from ref. |38) where n + d 
scattering with a simple model potential has been evalu- 
ated numerically for energies up to 1500 MeV and com- 
pared with the Glauber approach. The path-integral rep- 
resentation and this variational approach will certainly 
be also useful in attempts to evaluate the scattering am- 
plitude by stochastic (Monte-Carlo) methods [39^,10 . Fi- 
nally, as the present approach relies solely on actions and 
not wave function(al)s rendering it a priori more general, 
it seems that applications to nonperturbative field theory 
may also be possible. 
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6 (degrees) 

Fig. 3. (Color online) Same as in fig. [T]but for the various approximations by Chen and Wallace's 3rd order eikonal expansion. 
For comparison also the best variational result (ABc 31 in the naming scheme (14. lip ) is shown. 
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Appendix 



A. Evaluation of averages for the variational 
calculation 

In this section, we evaluate the different averages entering 
the Feynman- Jensen variational principle. 

We will be using a short-hand notation which simpli- 
fies the algebra considerably: first, as usual the convention 
that over repeated indices is summed and second that con- 
tinuous variables (times) are treated the same way. More 
specifically, i,j ■■■ = 1,2,3 denote the cartesian compo- 
nents of vectors and matrices, a, /3 . . . = 1, 2 ... (3 -I- d) 
the components of super-vectors and A = {a,t) etc. com- 
bine discrete and continuous indices. V^,/^'^... denote 
the transpose of the corresponding quantity. We do not 
distinguish between upper and lower indices. Occasionally 
we also use the dot to denote a scalar product between vec- 
tors or between super-vectors. Thus, for example the trial 
action (13. 7p may be written in several different forms as 



/-t-oo u-ru 
dtdt' J2 ^<.it)Ao,f3it,t')Vpit') 
a, = 1 

/+00 3+d 
dt Bait) Yc.it) 

- OO „ 1 



= -V„(i) A„^(i, t') ¥p{t') + B„(t) V„(t) 

TTl 

= —V ■ AV + B ■ V 



(A.l) 



showing the obvious advantage of a concise notation. 
We first need the path integral over the trial action 



mo 



J V^vV^w exp ( iAt ) (A.2) 



which from the rules of Gaussian integration can be im- 
mediately computed as 



mo 



constant 
Be^/^iK) 



exp 



2m 



(A.3) 



Here 



DetX = exp[Tr \nX] 



(A.4) 

denotes a determinant, both in continuous (also called a 
functional determinant) and discrete variables and con- 
sequently "Tr" is a trace in all variables. The path in- 
tegral in eq. (I2.4ap is normalized such that it is unity 
without interaction, i.e. B = and A = 0-3. Here 0-3 is 
the (extended) third Pauli matrix which appears because 
the kinetic term of the anti- velocity has an opposite sign. 



i.e. (T3 = diag(l, 1, 1, — 1, — 1, — 1) for the "aikonal" rep- 
resentation which uses a 3-dimensional anti-velocity and 
(T3 = diag(l, 1, 1, — 1) for the case of a 1-dimensional w. 
This normalization determines the constant as Det ' 0-3. 
Defining for convenience the vector 



I.e. (La 

f +00 



or 



Ca{t) 



dt' K~lit,t')V>fiit') (A.5) 



and using the fact that is its own inverse we finally 
have 



mo 



= Det"^/2(cr3A) exp ( - 



exp 



2m 



-iTrln(a3A) - -^d 
2 2m 



(A.6) 



This will also serve as a master integral or generating func- 
tion for the averages we have to calculate. 

The Feynman- Jensen variational principle requires the 
average of the difference between the full action and the 
trial action 



{AA), ^ -(V.(a3-A)V),-(] 



(A.7) 



Here 

x(b,v] 



/+00 
dtVi^.^fit) + 
Equant it)) (A.. 
-00 



is the phase in the path integral (I2.4ap for the impact- 
parameter 5-matrix and the average has been defined in 
eq. ([331) ■ 



A.l Computation of (B V)^ 

The evaluation of this average is easily done by putting 
an artificial scalar factor a in the linear term of the trial 
action and differentiating with respect to this factor 



(B- V)^ 



-i — In mo(C 
da 



2m 



da 
m 



(A.9) 
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A.2 Computation of f ( V • (ag - A)V)t 

We obtain this average by the same technique: differentiating mo with respect to B brings down a factor V. Therefore 



TO 



-(V.(a3-A)V), 



1 m S^niQ 



After the insertion of 



J^TOo 



Too 



I 

TO 



we obtain 



y(V-(a3 



or, in continuous notation, 



TO 

T 



TO 



y (V-(a3-A)V), = ^[C-a3C-C-AC] + ^Tr[^3A-i-l] . 
Together with eqs. (jA.6p and (jA.9p we thus have for the variational functional 

S(b,A,C) = expi -Lc.a3C-iTr[ln(a3A)+a3A-i-l] +i (x)^ 
2m 2 



(A.IO) 



(A.ll) 



(A.12) 



(A.13) 



(A.14) 



A. 3 Computation of {x)t 

Finally the average over the phase (lA.Sp can be done by a Fourier transformation of the potential: 



+ 00 



dt 



(27r) 



V^(p) exp [ip • Xref(i)] (exp (zp • Xquant(0 ) )t 



To simplify the calculation we write 

Xqua„t(i) =: rit,t')Yit') ,i.e. x^,,„t(t) 



dt' r,a{t,t')Ya{t') 



(A.15) 



(A.16) 



where riii{t, t') is a 3 x (3 + (i)-matrix operator. Note that as a messenger between ordinary space and the "super" -space 
in which velocity -I- anti-velocity live it cannot be a square matrix. Indeed, from the connection (|2.6p between quantum 
fluctuation and velocity we immediately find 



. /sgn(i-i') 
r{t,t') = -\ sgn(i-0 



-sgn(-t') 
-sgn(-t') 



sgn(i - t') 



-sgn(-t')> 



(A.17) 



for the "aikonal" case and from cq. (j2.9l 



-, /sgn(i-<')-sgn(-t') 

r{t,t')^-\ sgn(i - i') - sgn(-t') 











sgn(t-t') -sgn(-t'), 



(A.18) 



for the "ray" case if we choose the mean momentum along the 3-axis. 

In any case the average in eq. (jA.151) involves only Gaussian integrals where in the exponent the 



term linear in V = 



i [Ba + PiFia ] Va for the numerator . 



iB • V 



for the denominator. 



(A.19) 
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Thus 



;exp(ip-Xq„ant(i)) )t = ^ J T^^'T^'^w exp|^^V^AF + ^[B + p•^(^,•)] • v| 



2m 

xexp<j ^B-A^^B \ , (A.20) 

where the last factor comes from the normalization of the average by l/mp and no summation/integration over the 
external time-parameter t is implied. Working out the exponential we find terms linear and quadratic in p while the 
constant term cancels. The term linear in p gives the trajectory 

pS) xJrf(i)-^r,„(t,t')A-i(t',0%(^"), (A.21) 



or in short-hand notation 



The term quadratic in p reads 



p{t) = x,ef(i)--r(t,-)C. (A.22) 
m 



{t, t') A-i [f, t") r,p (t, t") p, (A.23) 

which becomes —{l/2)piPjaij(t) if a square (3 x 3)-matrix 



<j,j{t) = -r,^{t,t')r,0{t,t")A-^^{t',t"), or <j{t) ^ - r{t,-)A-^r^{-,t) (A.24) 



is defined (again no summation/integration over the external t). From the definition p.l4p for the Gaussian transform 
of the potential we then obtain the following simple expression for {x)t ■ 

r + oo 

{X)t^~ dtV^itMt)). (A.25) 



B. Variational equations 

In app.|S]we have calculated the quantities to be varied in the Feynman- Jensen functional. The result (jA.14p may be 
written as 

5(b, A, C) = exp [ iXo + iXi - Q{(j^k) ] , (B.l) 

with 

/ + 00 

dtV,(t){p{i)), (B.2) 
-oo 

Xi := -^C • fiaC, (B.3) 
Zra 

Q{X) ■.= ]^Tr[\x\X + X-^ -l] . (B.4) 

The argument of the potential - the trajectory - is given in eq. (jA.22[) and its Gaussian width in eq. (jA.24[) . We will 
now compute the variational equations for A and B which follow when stationarity of ^(b, A, B) is required. However, 
since only the combination C = A^^B (but not B) and A~^ enter the functional to be varied, it is more convenient to 
vary with respect to C and A~^. 
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B.l Variational equation for C 

The vector C only shows up in the phase Xi and the trajectory p{t). Therefore stationarity of the variational functional 
requires 

^-£..a.K,.,WO)|§io ,B,5, 

where we have used the chain rule in the last term. From eq. (jB.Sp we immediately have 

6Xi 1 



and from eq. (|A.22p 



('T3C)„ (t) (B.6) 



= — r,a{t\t). B.7 

d'La(t) m 

Thus 

(a3C)„(0 = ^d,V,{t')r,^{t',t) . (B.8) 
As (73 is its own inverse and diagonal in the time variables, we may solve this as 

C„W = -{<Ji)^pd,V,{t')r,p{t\t), (B.9) 

or 

Cvar = -fTaT^J. (B.IO) 

where 

J.W = ^v^itMt)) (B.ii) 

is the vector of potential derivatives (also called the Jacobian). We can now put this expression into the trajectory to 
get 

p,,,(t) = x,rfW ~ lr(t,s)C™'^(s) = XrefW + -r{t,s)a^r'^{s,t')3^{t') (B.12) 
m m 

=:x,rf(t) + i:o(i,OJ^(i')- (B.13) 
Here we have defined the (3 x 3)-matrix 

Z'o(i,t') -r{t,s)a^r^{s,t') (B.14) 
m 

whose explicit form can be worked out with the help of eqs. (|A.17I) and (lA.lSp and the relation 

/ ds sgn{s-t)sgTL{s~t') = 2[T -\t-t'\] . (B.15) 

J-T 

Due to the contribution from the anti-velocity all terms which would diverge for T — > cxd cancel and one obtains the 
expressions given in eqs. p.lSp and p.l9ap . Note that is the Green function for the corresponding equation of 
motion. For example, we have 

so that one can write 



^ (4^-^)) Jt,0 = '5(t-t')<5., (B.16) 



4'"'^ = ^ ^Jfl (B.17) 
—mof 



where V denotes the principal value prescription. Indeed, evaluating in Fourier (E-) space with {E\t) = e:xjp{iEt) / \/2tt 
as transformation function we obtain 

exp{iE{t ^ t')) _ dij 1 /■+°° ,^ exp(i£;(< -t')) - 1 
/ dE 



in agreement with eq. p.lSp . 





+ 00 








-oo 






dE 


TO TT Jo 





E^ m 271 J_^. E^ 

1 - cos(E{t - t')) 1 



E^ 2m 



\t~t'\S^j, (B.18) 



J. Carron, R. Rosenfelder: Variational Approximations in a Path-Integral Description of Potential Scattering 



19 



B.2 Variational equation for A ^ 

The matrix A^^ only enters in the Gaussian width a (appearing in Xq) and in the "fluctuation term" Q. Therefore 
requiring stationarity means 



5n . /■+°°^^ 5V„{t) 



From the definition p.l4p we find 



dt — ' ^ = . (B.19) 



= ^d.d,v^{t) ' . (B.20) 



Thus, we have to compute the functional derivative of the Gaussian width and fluctuation term with respect to A ^. 
This is easily done for the Gaussian width: From eq. (IA.24[) it follows: 

^""^'^^^ - -rutM)rjp[tM) ■ (B.21) 



For the fluctuation term we rewrite eq. (jB.4p as 



Q{a^k) = ^Tr [ - ln(A- V3) + (A" V3) - 1 ] (B.22) 



and use the chain rule (in our condensed notation from app. 



1 - ;)rA-i^„i„„ — - -WI^YZT— Sac {<J3)bd ■ (B.23) 



^^^g 9[A-i(73]cD aA^B d[k-^a3]cD 

Since 

dTvfiX) 



„ = [f'iMDC (B.24) 



^CD 

for any differentiable function f{X) we obtain 

^^ll^ = -\{[[^-'<yA'') ^^-5ad} 5Ac[cr,)^^ = -]^{kAB-[<y,]AB), (B.25) 
where we used the fact both and A are symmetric. In extended notation the resuhs reads 

fl , = -\{^c.p{tiM)-{^,)^,5{ti-t2)] ■ (B.26) 
Utilizing eqs. (|B.20p . (|B.21|) the variational equation (jB.19|) therefore becomes 

1 r 1 1 r^°° 

-^[ka,0{h,t2)-ia3)^pd{h-h)j = — J dtd,djV^{t)r,^{tM)rj0{t,t2)- (B.27) 
In matrix notation this simply reads 

Avar = a3--r^H„r (B.28) 

TO 

where He stands for the matrix of second derivatives (also called the Hessian) of the potential 

WJ{t,t') :- d,dj K.(t)(p(t)) S {t - t') . (B.29) 
It is possible to put the results in an even more compact form by defining 

S := —rk-^r^ . (B.30) 

TO 



This can be easily proved by a power series expansion of f{X), piecewise differentiation and resumming. 
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Note that for the free case A = 0-3 this reduces to the matrix Uq defined in eq. (IB.14P which is the Green function for 
the equation of motion (jB.lSp . Multiplying eq. (IB.28P from left by and from right by A^^F'^ /m, the variational 
equation for then translates into a Lippnrann-Schwingcr-like integral equation 



(B.31) 



All relevant quantities can be written in terms of S{t,t') and the variational trajectory p{t). For instance from 
eqs. (|B.3p . (IB.lOp and (jB.14p we immediately obtain 



Xi — 2 '^'^ ' -^0 Jo 



(B.32) 



Less obvious is that the fluctuation term fi^a^A) defined in eq. (jB.4l) may also be expressed in terms of U and p. 
However, if the variational solution (jB.28p is multiplied by and inserted into 



I2(ct3A) = --Tr [ln(A"V3) - A^Vs + 1] 



(B.33) 



one obtains 



o — Tr 



In 1 + A 



-1 1 



(B.34) 



Using the definition (IB.30p and the properties of the trace this becomes 



n = --Tt[ ln(l + SH,) - SH, ] 



(B.35) 



Finally from eq. (jA.24p we see that the Gaussian width a{t) is just given by the diagonal element of lE. 



C. Calculation of the second cumulant 



In this section we compute the first correction to the variational approximation, i.e. the second cumulant 



[AAf) -{AA) 



AA= -Y- ((73 - A) ' 



(G.l) 
(C.2) 



for both representations. Since the computational steps are similar to the ones in app. El we will not enter into too 
much details. After some algebra we find that the second cumulant is given by the following terms which we will 
evaluate in the subsequently: 



A2 = (x'),-(x)? + -C-AC + z,n-f 



m Stuq S{x)t , 
2 — r^-::^ — rrr + m- 



a=l 



{Ea)t + -{Xa)t 

m 



mo (5Ba(t) 6Mfs{t') 
where the quantity E stands for 



+ 



SHE) 



,{t)5Mp{t') 



{E), := (V.((T3-A)V), . 
From our study of the first cumulant, eqs. (|A.13p and (|A.25p . we can compute 

, , 1 



da 



V Xa , 



• rc. 



and 



d_ 

da 



a=l 



(C.3) 



(C.4) 



(C.5) 



(C.6) 
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Similarly, we have 



1 (5mn 



mo SMa{t) 
Hx)t 



= C„(0 

m 



-\wVa-rA-^] (t) 



HE), 



while 



(5B„(t)<5B^(t') 



Q/3 



1 

m 
1 



2 [A-'r^H,rA-']^^it,t') 



(C.7) 
(C.8) 
(C.9) 

(C.IO) 
(C.ll) 
(C.12) 



Here, the "var" above the equality sign in the last line means that the variational solution (IB.28P has been used. We 
have now everything needed to compute the second cumulant. After some (tedious) algebra we obtain 



A2 = (x'),-(x)?--C-a3A-V3C + iTr 



^3 - 1] - Tr [A-V3 - 1] . 



Using the variational equations and the Lippmann-Schwinger-like equation we find 

Tr[(73A-i -1] =Tr[SH,] 



-Tr[(73A-V3A-i - 1] =Tr 



while 



so that our final result is 



— IL • (T3A (T3II 
m 



(C.13) 

(C.14) 
(C.15) 

(C.16) 
(C.17) 



The quantity ( ) j is evaluated in the very same way as ( X )t j i-^- one transforms the potential to Fourier space, and 
performs the functional integrations. The result is the one given in eq. p.27p . 

D. Some numerical details 

Here we discuss several points which are essential to obtain reliable numerical values for our variational approximations. 



D.l Variables 

As the range of the potential greatly determines the dynamics of the scattering process we switch from time variables 
to distances by substituting 

t = ^ (D.l) 

^char 

where the characteristic velocity is chosen as the asymptotic velocity on the different reference trajectories, viz. 



Vcha 



K_ 



"aikonal" 
"ray" . 



(D.2) 



This is reasonable for high-energy scattering where indeed the particle mostly travels along the reference trajectory 
and implies that each power of t and each integration over t is suppressed by an inverse power of K (or k). It becomes 
less convincing at low energy where the characteristic velocity may be totally different from the asymptotic one - for 
example, in scattering via a resonance. However, all our applications will be in the high-energy domain. 
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D.2 Numerical integration 

Initially our numerical integrations have been performed 
by using Gauss-Legendre quadrature rules for a finite in- 
terval [— Zmax, +-Zmax] • Typically we have taken 



(4^6) -i? 



(D.3) 



as we expect that outside this interval the potential is 
practically zero. Of course, Zmax had to be varied to ensure 
stable results. In previous work we had mapped the infi- 
nite z-interval to a finite one (for example by z = i? tan ip) 
but in the present case this caused some problems for the 
iterative solution of the variational equations which are 
avoided if the integration is over a finite, not too large 
interval. 

The numerical integration was performed with n^. sub- 
divisions and Ug Gaussian points in each subinterval re- 
quiring N = Tie X Hg function calls altogether. Again 
the number n,, was varied to verify stability of the re- 
sults whereas the number rig was kept fixed at moderate 
values which is advisable for oscillatory integrands, e.g. 
(ugjUe) — (32,2). Since we deal with multi-dimensional 
integrals and have to solve the variational equations for 
each b-value on the grid the computing time rapidly in- 
creases when making the grid finer and finer and at large 
scattering angles (where the scattering amplitude becomes 
small by interference) it became difficult to obtain sta- 
ble numerical results. In these cases use of the adaptive 
routine DCUHRE [ITll^ for the b-integration was quite 
helpful as it chooses the integration points according to 
their relative importance. 

Later we realized that integrating numerically over 
nonanalytic functions like \z—z'\ or \z'\ (which occur in the 
Green function So{z, z') ) is not well handled by Gaussian 
or similar quadrature rules. The reason is that their error 
is proportional to some high derivative of the integrand 
within the integration interval which makes them suitable 
for analytic functions but not for nonanalytic ones. Conse- 
quently, the simple trapezoidal rule (see, e.g. eq. (25.4.2) 
in ref. [13], with h = {zn — zq)/N denoting the increment) 



dz f{z) — h 



I + A + 



+ fN-1 + -J- 



12 



no , zo<(<ZN (D.4) 



is as efficient (or better: inefficient) in integrating nonana- 
lytic functions as a A^-point Gaussian integration which is 
exact for polynomials up to z"^^^^ or Simpson's rule whose 
error is proportional to h^f^^\0- This is demonstrated in 
table m where the test cases 



+ 00 



h{z): 



dz' \z — z'\e 



\pKZ ■ erf (z) 



•(D.5) 



+ 00 



h{x) y dz' [\z~- z'\-\z\-\z'\^e-''" 

— OO 

= /i(z)-V^|z|-l (D.6) 



are considered as typical examples of integrals for a Gaus- 
sian potential in the "aikonal" and the "ray" represen- 
tation, respectively. Here erf(z) denotes the standard er- 
ror function with erf(— z) — — erf(z). It is obvious that 
the accuracy is relatively poor and the convergence with 
increasing number N of function calls is disappointingly 
slowEI. This can be improved in the following way: 

For simplicity, we first consider integrands of the /i-type, 

f{z') = \z-z'\cl^{z')+g{z') (D.7) 

where (f>{z'), g(z') are differentiable functions which van- 
ish rapidly enough for z' ±oo. Therefore to a good 
approximation 



dz 



'{|z-z'|</)(z') + .g(z')} 



dz' 



(z - z') <^(z') 



+ 2„ 



dz' 



{z'-z) 4>{z')+g{z') . (D 



Note that the integrand is analytic but different in the two 
integrals. We may now apply the Euler-MacLaurin sum- 
mation formula (ref. [43 , eq. (25. 4. 7)) assuming that the 
point of nonanalyticity z — kh in a, multiple of the incre- 
ment h = 2zniax/-A- As we are evaluating the variational 
equations iteratively on a grid this certainly is the case 
and we obtain 



dz 



(z') -I- g{z') I ~ trapezoidal rule 



E 



2m I J2m — 1 



(2m)! 1 dz'2"-i 



(z-z')0(z') + .g(z') 



dz'2'»-l 



(z'-z)0(z')+5(^') 



/l2 ^4 

trapezoidal rule + —cpk - -t:;:^!' + O (h^) . (D.9) 
6 120 

where i?2m are the Bernoulli numbers. Note that the cor- 
rection terms come from discontinuities of derivatives of 
the integrands below and above the point z — kh (we 
neglect the contributions around the endpoints at iz^ax 
since the integrand is assumed to be very small there). 
The integration rule up to 0{h^) is denoted by "Euler- 
MacLaurin2" in table |H Keeping correction terms up to 
0{h'^) ( "Euler-MacLaurin4" ) requires the second deriva- 
tive of (/)(z = kh) which we simply approximate by 



1 

ft2 



'Jk-1 



- 20fe 



(D.IO) 



For comparison: without absolute sign in the integrand, i.e. 
for an analytic function, the absolute deviation from the exact 
result is lower than 10~^^ already at A'^ = 48 for all integra- 
tion rules. Strictly speaking, of course, the equally spaced rules 
require + 1 function calls. 
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Table 4. Maximal deviation of the numerically evaluated integrals Ii/2{z) (see eqs. (|D.5|) . (|D.6|) ') from the exact value for 
different quadrature rules and function calls A'' in the interval z G [—5,5]. The integration range was made finite by cutting 
off the integrand for \z'\ > 5. "Gauss 24" etc. refers to a Gauss-Legendre quadrature rule with Ug — 24 points etc. The last 
two lines give the results obtained with the Euler-MacLaurin summation formula retaining corrections up to second and fourth 
power in the increment, respectively, due to the nonanalytic behaviour of the integrand. 



Ii{z) 



Hz) 



Rule N ->■ 

i 


24 




48 




72 


144 


24 


48 




72 


144 


Gauss 24 


6.38 ( 


-2) 


5.12 ( 


-3) 


7.53 (-3) 


1.24 (-3) 


1.01 (-1) 


5.12 


(-3) 


1.14 (-2) 


1.24 (-3) 


Gauss 48 






1.72 ( 


-2) 




1.94 (-3) 




2.61 


(-2) 




2.91 (-3) 


Gauss 72 










7.77 (-3) 


5.80 (-4) 








1.17 (-2) 


5.78 (-4) 


Simpson 


5.04 ( 


-2) 


1.40 ( 


-2) 


6.34 (-3) 


1.60 (-3) 


5.30 (-2) 


1.42 


(-2) 


6.37 (-3) 


1.60 (-3) 


Trapez 


2.95 ( 


-2) 


7.27 ( 


-3) 


3.22 (-3) 


8.04 (-4) 


-2.95 (-2) 


-7.27 


(-3) 


-3.22 (-3) 


-8.04 (-4) 


Euler-MacLaurin2 


5.25 ( 


-4) 


3.17 ( 


-5) 


6.23 (-6) 


3.88 (-7) 


-7.57 (-4) 


-4.59 


(-5) 


-9.01 (-6) 


-5.61 (-7) 


Euler-MacLaurin4 


6.33 ( 


-5) 


1.00 ( 


-6) 


8.81 (-8) 


1.39 (-9) 


-9.91 (-5) 


-1.60 


(-6) 


-1.43 (-7) 


-2.23 (-9) 



Table 2] shows that this gives vastly improved numerical 
results for the test functions: two or three orders of mag- 
nitude more accurate than the simple trapezoidal rule. If 
there is another point at z' = (as in the Green func- 
tion of the "ray" representation or in the integrand of 
hiz)) we simply add the corresponding correction for that 
point. For the variational calculation we used the "Euler- 
MacLaurin2" integration rule with = 60 — 120 integra- 
tion points. 



or e = 10~^ which results in something like a dozen itera- 
tions for small b whereas only one iteration is needed for 
large (peripheral) values of the impact parameter because 
the potential for those trajectories is already very feeble. 
The criterion (ID.12[) makes sense because the integrand is 
proportional to ex-p{i{Xo+Xi -f ...)) — 1 so that consider- 
able computing time is saved for trajectories which barely 
feel the influence of the potential. 



D.3 Variational equations 

We have solved the variational equations by iteration start- 
ing with the free solution. Updating during iteration is 
performed by the simple "linear mixing scheme" 



'^mix ?/out + (1 ~ \nix) U; 



in) 



(D.ll) 



although more elaborate schemes are available (see e.g. ref. 
[33] and references therein) . Simply taking Amix — 1 works 
quite well for high-energy, forward scattering since the po- 
tential contribution is suppressed by appropriate powers 
oi K = fc cos (9/2 in the "aikonal" representation, or k in 
the "ray" representation. In the former case convergence 
obviously deteriorates at backward angles and one needs 
more and more iterations to fulfill the requirement that 



(^0+^1 



sin) 



< e 



(D.12) 



which is imposed at fixed impact parameter b before the 
iteration is allowed to stop. Typically, we take e — 10~^ 



D.4 Gaussian integrals and transforms 

Although in all textbooks the Gaussian integral in n di- 
mensions is given as 

In iy) = J rf"a; exp [-x'^Cx + iy'^x] 



t"/2 



Vdet„ C 



exp 



, (D.13) 



this strictly holds only if the matrix C is real symmetric 
(or hermitean) positive definite. Recall that this property 
implies positive eigenvalues Xj so that the determinant 
expressed as the product of the eigenvalues 



det„C = Jl > 

3 = 1 



(D.14) 



is positive and eq. (|D.13I) is unambigous. However, in our 
case 

C = A + iB (D.15) 
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is only complex symmetric (the quadratic form x^Cx 
projects out the symmetric part) and for convergence of 
the multidimensional Gaussian integral the real part A 
has to be positive (semi-)definite, i.e. x'^Ax > for all 
X. Under these conditions eq. (ID.lSp still holds provided 
the correct sign of the complex square root is chosen. This 
is a subtle but important point about which we haven't 
found very much in the literature - except vague remarks 
that "the sign of the square root is fixed by ... analytic 
continuation" (ref. [45], p. 421). A possible sign change 
is equivalent to an additional phase tt in the exponent of 
eq. (jP.lSp and thus very similar to the Maslov phase cor- 
rection (multiples of 7r/2) for a semi-classical propagator 
when the trajectory of the particle goes through a focal 
point pS] . 

In the following we outline how a proper treatment 
of these "branch corrections" may be obtained. First, we 
specify which complex square root (which has branch- 
points at and oo) we will use in the following: we choose 
the cut between these branchpoints along the negative real 
axis and define the principal square root of a complex 
number z as the one with a positive real part, i.e. 



words: we can write 

Xj = Tj e"^^ , with < 



'z := \/ \z\ exp [ i i ar| 



argz < TT 



(D.16) 



This is also the value returned by the subroutine for the 
complex square root in the numerical implementation. Then 
we assume that the complex symmetric matrix C may be 
diagonalized by a complex orthonormal transformation 



C = ODO^ 



with O^O = OO^ = 1 (D.17a) 



and 



so that 



D = diag(Ai, A2 . . . A„) 



(D.17b) 



Iniy) = J d"x' exp 



n 



dx" exp 



^Xjx'- - 



4A,- 



(D.18) 



because the Jacobian of an orthonormal transformation is 
unitv F°l. Due to a theorem by Bendixson (see eq. (6.9.15) 
in ref. [17]) the real parts of the complex eigenvalues 



ReXj > 



(D.19) 



remain nonnegative when the eigenvalues of the real part 
of C are assumed to be (semi-)positive definite. In other 



Lacking an explicit proof we here assume - as done tacitly 
in all textbooks - that the complex orthonormal transforma- 
tion of the original co-ordinates Xi and the linear shift leads 
to integration paths for the transformed co-ordinates x" in the 
complex plane which can be safely rotated back to the real axis. 
Analytic continuation of the real result faces a similar difficulty 
as one can only perform it from a region in the complex plane 
and not from the real axis. 



(D.20) 



Thus, each integral in the product of eq. (jP.lSp is conver- 
gent and we obtain 



exp 



1/2 



1 



exp 



1 " 



{\{^ exp(-iy^C7-iy) (D.21) 

as expected. However, it is important to note that the 
prefactor is not given by the inverse square root of the 
determinant of C but by the product of the inverse square 
roots of the complex eigenvalues (ID.20p . This may be dif- 
ferent [3 depending on the value of 



IL 



(D.22) 



Then 



n 



\ yXj ^/AetnTJ 



exp ( —i nbr(C') tt ) 



,, ^ I 6xp ( -■^argdet„C-inbr(C)7r )(D.23) 
V|det„C| V ^ 



where the "branch number" is given by 

|<Z>(") I -f TT ■ 



nbr(C) 



27r 



(D.24) 



and [x\ is the maximum integer not greater than x. 
"■br = 0,2... denotes the principal branch of the square 
root and ribr = 1,3... the other, negative branch. From 
eq. (ID.22P we have the following bound for a. {nxn) com- 
plex symmetric matrix C with nonnegative real part 



nbr(C) < 



(D.25) 



Note that in contrast to the Maslov correction this ad- 
ditional phase is not discontinous but rather corrects the 
phase jump when crossing the (arbitrary) branch cut of 
the complex square root. In this way the analytic continu- 
ation of the result (ID.13|) is achieved. As an example take 
^ = TT + e to obtain — ^arg det C — ribr tt — — 7r/2 + |e|/2 
for e < (ribr — 0, above the cut) and tt /2 ~ e/2 — tt = 
— 7r/2 — e/2 for e > (ribr — 1, below the cut). 



Take the simple example: n = 3, Aj = A = i -f 0^. Then 
deta = -i-0+, ^^^iet^ = (l-i)/^but (v^A)^ = -(1-i)/^. 
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After these preliminaries we can now calculate the gen- 
eral Gaussian transform required in our variational ap- 
proach. This is particularly straight-forward and simple 
for a Gaussian potential: from 



K(p) - V{p) -exp 



1 



exp 



2a,y /i?2) 



(D.26) 



it follows by inverse Fourier transformation 



K(x) = Vb^'/'i?' 

Vr 



exp 



Pi ^ij Pj ^Pi^i 



(27r)3 

° exp [—aXiC^j^ Xj — mbr(C)7r] (D.27) 



where 



</det3 C 



2a a 



1 



(D.28) 



is a complex symmetric (3 x 3)-matrix whose real part 
should be positive (semi-)definite. Unfortunately we were 
unable to veriiy this property analytically for our varia- 
tional solutions but did not encounter any numerical in- 
stabilities (which would be caused by a blow-up of V^) 
during the iterative solution of the variational equations 
if -s^max was not too large. Determinant and inverse of the 
(3 X 3)-matrix C are known from elementary calculus but 
we have not found a simple but reliable method to de- 
termine the "branch number" nbr which in principle - 
according to eq. (jD.25|) - could be nonzero even in this 
case. A method of "branch tracking" has been described 
in ref. [15] but we used a less elegant, brute-force approach 
in which the complex eigenvalues were determined numer- 
ically with the NAG routine F02GBF and the prefactor 
was calculated as a product of the square roots of these. 
As a by-product it was confirmed that the real parts of 
the eigenvalues were always positive. After evaluation of 
the Gaussian transform of the potential the Jacobian and 
the Hessian then simply follow by differentiation: 

(J,),, = d, K(x) - -2a (C-i)^^ Xk K(x) (D.29) 



didj Vct(x) 



-2a 



(D.30) 



-2a (C-i)^^Xfe (C-i) K(x).(D.31) 



For the calculation of the second cumulant we also need 
the double Gaussian transform 

/g := J (fpid^p2 exp[~p^Cp + ix-p] (D.32) 



where 



C 



Cii Ci 



(D.33) 



Ci2 C22 

is a complex symmetric (6 x 6)-matrix (see eq. p.27p ) and 

(D.34) 





, X = 











are 6-dimensional (column) vectors. Determinant, inverse 
and branching number have been evaluated as in the 3- 
dimensional case by calculating the complex eigenvalues 
- a procedure which increased the execution time of the 
program considerably. We found that (under our kine- 
matical conditions) no branch-crossing occured in the 3- 
dimensional calculation of the Gaussian transform of the 
potential but the sign-change of the complex square root 
was essential to obtain the correct results in the 6-dimen- 
sional case (calculation of the second cumulant). 



D.5 Calculation of Q 

Arising from a functional determinant the quantity i? 
poses a particular problem for numerical evaluation. Sev- 
eral approaches are possible: 

First, one may employ the classic method of Gel'fand 
and Yaglom (already contained in textbooks, e.g. in 
ref. [33], Chapt. 6) to calculate a functional determinant 
as solution of an initial value differential equation. Indeed, 
defining 

^0 ^Tr \n{l-S^H„) (D.35) 

(so that Q = Qq + Ti-{S Her)/ 2 ) its exponential is given 
as ratio of two functional determinants 



exp (2 Qo) 



Det {-df - H^/m) 
Dct(-a2) 



(D.36) 



and one may apply the Gel'fand- Yaglom procedure to 
evaluate it. However, the boundary conditions for the eigen- 
functions f{t) of the differential operators in eq. (jD.36l) are 
not of Dirichlet type as in the standard method but (in 
the "aikonal" case, cf. eqs. (|3.32ap and (I3.32b[) ) of the 
form 



lim 

T-J-oo 



= 



{/(r) + /(-T)-r(/(t)-/(-r))} 

{/W + /(-r)} =0.(D.37) 



lim 



Although Kirsten and McKane [5^ recently have general- 
ized the classic procedure to more general boundary con- 
ditions hke (|D37)) we do not follow this approach since 
simpler alternatives are available. 

These include, second, the calculation of f2 as power series 
in either SH^ or S^H„: 



■Tr 



2 

1 °° 
iTrV 

2 ^ 

n=2 



In (1 + SH, 

i-r 



Tr {SH„ 



[2 = 



iTr 
2 

1 



-TrV 
2 ^ 

n=2 



1 



1 



Tr (Z-oiJ,.)" . 



- 1 



(D.38) 
(D.39) 

(D.40) 

(D.41) 
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Table 5. Comparison of values for the fluctuation term (functional determinant) fi calculated in different ways (see text) at 
b/R — 1 for the Gaussian potential (|4.1|l . Parameters of the calculation: Zmaic/R = 5, (ng,ne) = (32,2), e = 10~^, Amix ~ 1- 
Results for two values of the energy and the potential strength are displayed to allow comparison with the "aikonal" phases to 
which the total sum Xo + Xi -\- ifi should tend in the high-energy limit. 



KR = 4, 2mVoR^ = - 
( # of iterations = 5 



Re 



Im 



KR = 8, 2mVoR^ = -8 
( # of iterations = 4 ) 



Re 



Im 



n , eq. (|D.39|| 


-4.1931 (-3) 


-6.2753 (-4) 


-9.6244 (-4) 


-4.3797 (-5) 


n , eq. (|D.41|| 


-4.1930 (-3) 


-6.2751 (-4) 


-9.6258 (-4) 


-4.3810 (-5) 


n , eq. (|D.43|| 


-4.1935 (-3) 


-6.2748 (-4) 


-9.6243 (-4) 


-4.3812 (-5) 


n , eq. (|D.38|| 


-4.1934 (-3) 


-6.2747 (-4) 


-9.6246 (-4) 


-4.3814 (-5) 


Q , eq. (|D.42|| 


-4.1932 (-3) 


-6.2756 (-4) 


-9.6247 (-4) 


-4.3799 (-5) 


n , eq. (|D.44|| 


-4.1932 (-3) 


-6.2756 (-4) 


-9.6247 (-4) 


-4.3799 (-5) 




3.4056 (-1) 


8.1621 (-3) 


3.3382 (-1) 


1.8965 (-3) 


Xo+ Xi+ifl 


3.4118 (-1) 


3.9688 (-3) 


3.3386 (-1) 


9.3403 (-4) 


(0) 1 (1) 1 (2) 

xVi + xVi + Xai 


3.4176 (-1) 


0. 


3.3394 (-1) 


0. 


■ (2) 


0. 


3.3216 (-3) 


0. 


8.3041 (-4) 



If the variational equations are fulfilled this should give 
identical results as should the "mixed" form 



oo ^ 

V - Tr {SoHaT + Tr {SH,) 



(D.42) 



(D.43) 



in which 1/(1 — I^qHo-) — 1 = SH„ has been used. Note 
that in the "aikonal" case the sum also begins at n = 2 
since X!Q{t,t) = 0. If the variational equations are solved 
by iteration it is consistent to evaluate the various sums in 
eqs. (|D.39I) . (|D.41|) and (|D.43p up to n = # of iterations 
since each term is suppressed by an additional power of 
if or fc - if the iteration converges so will the power series 
expansion for f2. One should keep in mind that this pro- 
cedure only works as long as KR ^ 1 or kR ^ 1, i.e. at 
high energies and not too large scattering angles. 

Third, one may evaluate the functional determinant as 
a (3 X rig X rig)- dimensional ordinary determinant on the 
grid as one has discretized the time (or z-co-ordinate) for 
the solution of the variational equations and the various 
integrals anyway. The NAG program F03ADF was used 
for this purpose. 



Finally, the complex eigenvalues Xj of the matrix 1 -t- 
SHa^ may be calculated (with the help of the NAG routine 
F02GBF) so that from eq. (Id381) 



(D.44) 



This also checks whether a "branch crossing" may have oc- 
curred which is unlikely under these kinematic conditions 
as SH^ remains "small" and the power series of the loga- 
rithm is well converging. Indeed, we found no case where 
the square root of the determinant was different from the 
product of the square roots of the eigenvalues. 

Table [5] compares the results of the different methods in 
the "aikonal" representation at a fixed value of the impact 
parameter. One observes excellent agreement between the 
different methods. 

Also included is a test at high energies where according 
to eqs. p.46p and p.44p the phases and the fluctuation 
term can be described by the (much simpler) "aikonal" 
phases. These have been worked out in ref. ^22) for spheri- 
cally symmetric potentials (see eqs. (4.98), (4.99), (4.104) 
and (4.111) therein) so that the corresponding expressions 
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Table 6. Test of the Feynman-Hellmann relation (|D.49|) for fixed b/R = 1 in the "aikonal" representation and fixed bx/R ~ 
0.6, by/R = 0.8 in the "ray" representation. In the latter case the scattering angle has been fixed at 6 = 60° with the momentum 
transfer along the a;-axis. Accuracy parameters for the numerical solution of the variational equations are as in table [S] and the 
integration over the potential strength was performed by Gauss-Legendre integration with ufh points. 



"aikonal" 
Re Im 



"ray- 



Re 



Im 



r.h.s. of eq. (^49)) ufh = 8 3.41182 (-1) 3.96843 (-3) 
r.h.s. of eq. (|D:491) ufh = 12 3.41182 (-1) 3.96839 (-3) 



Xo + Xi + iO 



3.41183 (-1) 3.96882 (-3) 



2.50889 (-1) 3.38986 (-2) 
2.50889 (-1) 3.38987 (-2) 



2.50889 



3.38993 (-2) 



for a Gaussian potential read read 



D.6 Test of the Feynman-Hellmann theorem 



(0) 

Xai 



(2) 

Xai 



KR sy 2 

1 



{KR)"^ 16 V 3 



(D.45) 

(1 - 42/) e-'^y (D.46) 
1 - f 12 + V37r) y 

(D.47) 



(i2 + 2tt/V3 



-3y 



,(2) 
'^AI 



[KRf 



(1 - 4y + 2y 



2\ ^-2y 



(D.48) 



with C = 2mVoR'^/{KR) and y = b^/R^. 

Good quantitative agreement is observed which be- 
comes better at higher energies as expected. Note that 
at KR = 8 the strength has also been changed to keep 
mVo/K — constant. According to eqs. p. 461) . (I3.49P the 

difference between Xq + Xi + i f2 and J2k=o Xai + ^^Ai 
then should decrease as Indeed, a closer look at the 

numerical values in table [5] shows that this difference de- 
creases by a factor 7.9 in the real part and a factor 6.3 in 
the imaginary part when doubling the energy which is in 
reasonable agreement with the expected factor (8/4)'^ — 8. 



As explained in sect. 13. 61 the variational impact-parameter 
5-matrix should fulfill additional relations since its ingre- 
dients are solutions of the variational equations. Here we 
will use the dependence (j3.55D on the coupling strength 
Vq to test our numerical solutions. Upon integrating we 
should have 



Xo + Xi+iH 



dXjXo 



(D.49) 



for fixed b (and scattering angle). 

The r.h.s. of this relation may be evaluated by simple 
Gauss-Legendre integration (which has the advantage of 
avoiding the point A = 0) over the "phase" Xq at dif- 
ferent strength of the potential. As the integrand is a 
smooth function of the potential strength very few Gaus- 
sian points are necessary to achieve a stable result which 
is in excellent agreement with the l.h.s. of eq. (jD.49l) . 
This is displayed in table [HI both for the "aikonal" and 
the "ray" representation and constitutes a rather strin- 
gent test that our numerical scheme for solving the varia- 
tional equations and for calculating the stationary values 
is correct. 



Note addded in proofs. The correct mathematical frame- 
work to derive the Gaussian integral with complex sym- 
metric matrices is the theory of "pencils" [5T] (called 
"Biischel" in the German nomenclature [52]) of quadratic 
forms which avoids the problems indicated in footnote ^'^ . 
Fortunately this approach leads to the same result and 
procedures as used in Appendix D 4. 
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